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Abstract 

The aim of this study was the assessment of genetic diversity and differentiation patterns of European beech ( Fagus 
sylvatica L.) populations within its natural range in Poland. These stands cover 5.2% of the forest area in Poland, and 
form the predominant forest tree communities throughout the Carpathians and Sudety Mountains, and the moraine 
landscape of the Pomeranian Lake District. Varying environmental conditions have resulted in a great number of 
ecotypes and populations which are characterized by various ecological requirements. Poland represents the north- 
eastern limit of the beech’s natural distribution. 

Genetic diversity and differentiation was assessed using allozyme gene markers employing 9 enzyme sys- 
tems: glutamate-oxaloacetate transaminase (GOT - EC 2. 6. 1.1), leucine aminopeptidase (LAP - EC 3.4.11.1), 
isocitrate dehydrogenase (IDH - EC 1.1.1.42), malate dehydrogenase (MDE1 - EC 1.1.1.37), menadione reductase 
(MNR - EC 1.6.99.2), phosphoglucomutase (PGM - EC 2.7.5. 1), phosphoglucose isomerase (PGI - EC 5. 3. 1.9), per- 
oxydase (PX - EC 1.11.17) and shikimate dehydrogenase (SKDH - EC 1.1.1.25). 

Beech populations occurring toward the northeast of the natural range were characterised by a decreasing the 
average number of alleles per locus and percentage of polymorphic loci. The highest genetic differentiation was 
found in the East Carpathians. The genetic differentiation levels of Polish beech populations are comparable to the 
other European countries. 
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Introduction 

The aim of the present study was the assessment of genet- 
ic diversity of European beech (Fagus sylvatica L.) with- 
in its natural range in Poland. Beech is one of the most 
important forest tree species in Poland, and beech stands 
account for 5.2% of the total forest area (Lasy Panstwowe 



2009). They constitute the most typical forest tree as- 
sociations found within the lower forest belt in the Car- 
pathians and Sudety Mountains in the south of Poland, 
and throughout the moraine landscape of the Pomeranian 
Lake District in the northern of Poland ( Wojterski 1990). 

The present genetic structure of beech populations 
in Europe was formed by many different factors, inclu- 
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dig anthropogenic ones. Important factors affecting the 
gene pool were the glacial period and the locations of 
few beech refugia, from which beech colonized its pre- 
sent range during the postglacial, as evidenced by fos- 
sil pollen analyses (Szafer 1935; Raiska- Jasiewiczowa 
1983; Hazier et al. 1997; Huntley and Birks 1983). Other 
studies using genetic markers revealed existence of only 
one effective refugium (Gomory et al. 1999), suggested 
as the source of beech forests in the major part of Eu- 
rope except southern Italy and the Balkans on the ba- 
sis of the genetic differentiation monitored by isozyme 
loci, which was very weak. This hypothesis was con- 
firmed later by Magri et al. (2006) using a combination 
of genetic markers (isozymes and chloroplast DNA) and 
paleobotanical analyses (fossil pollen and macrofossils). 
Investigations of beech variation in Poland performed 
with isoenzyme studies (Sulkowska 2002; Gomory et 
al. 2003) showed high genetic diversity, similar to other 
neighboring European populations. Slight decrease of 
average number of alleles per locus and lower level of 
differentiation towards the northern natural range limit 
generally confirm the migration path after the glacia- 
tion period. 

Varying environmental conditions resulted in 
a great number of ecotypes and populations reflecting 
various ecological requirements (Dzwonko 1990; Gier- 
tych 1990; Sulkowska et al. 2008). In Poland the beech 
attains its north-eastern distribution limit (Szafer and 
Pawlowski 1972; Boratynska and Boratynski 1990). The 
geographic range of beech is limited by climate conti- 
nentality, soil conditions, winter temperatures and air 
humidity (Slawinski 1947; Jedlinski 1922; Boratynska 



and Boratynski 1990). The growth of beech stands out- 
side the natural beech limit indicates that the species 
possesses potentially wider range (Brzeziecki 1995; 
Tarasiuk 1999). 

Methods 

In this study 26 populations were sampled, represent- 
ing the best preserved selected seed stands, as well as 
a number of other natural beech stands (Tab. 1). Many 
of the beech populations chosen were located in the 
Carpathian Mountains, where beech is the most impor- 
tant forest tree species. 

Twigs and dormant buds were collected from 
50 trees in each population. Proteins from buds and 
cambium were extracted using 0.1 M Tris-HCl buffer 
pH 7.0. Starch gel electrophoresis was applied to gen- 
otype trees at allozyme loci. Histochemical staining 
and interpretation of zymograms followed Thiebaut 
et al. (1982), Merzeau et al. (1989), and Muller-Starck 
and Starke (1993). Nine enzyme systems were ana- 
lysed: glutamate-oxaloacetate transaminase (GOT 

- EC 2. 6. 1.1 - Got-2 ), leucine aminopeptidase (LAP 

- EC 3.4.11.1 - Lap-1), isocitrate dehydrogenase (1DH 

- EC 1.1.1.42 - Idh-1), malate dehydrogenase (MDH 

- EC 1.1.1.37 - Mdh-1, Mdh-2, Mdh-3), menadione re- 
ductase (MNR- EC 1.6.99.2 -Mnr-1). phosphogluco- 
mutase (PGM - EC 2.7.5. 1 - Pgm-1). phosphoglucose 
isomerase (PGI - EC 5. 3. 1.9 - Pgi-2), peroxydase (PX 

- EC 1.11.17 - Px-1, Px-2 ) and shikimate dehydroge- 
nase (SKDH - EC 1.1.1.25 - Skdh-1). Allelic frequen- 



Tab. 1 . Geographic location of Polish beech populations studied 



No. 


Provenance 


Geographic region 


Altitude (m. a.s.l.) 


Geographic coordinates (N/E) 


1 


2 


3 


4 


5 


1 


Cisna 


South (Mountains) 


900 


49° 11’ 


22 ° 18’ 


2 


Moczarne 


South (Mountains) 


800 


49° 06’ 


22° 29’ 


3 


Caryna 


South (Mountains) 


950 


OO 

O 

O 

OS 

'xt* 


22° 46’ 


4 


Rozsypaniec 


South (Mountains) 


950 


49 ° 05’ 


22° 50’ 


5 


Nad Sanem 


South (Mountains) 


750 


49° 15’ 


22° 42’ 


6 


Hulskie 


South (Mountains) 


650 


49 ° \ 2 ’ 


22° 38’ 


7 


Wielka Rawka 


South (Mountains) 


950 


49° 05’ 


22° 30’ 


8 


tosie 1 


South (Mountains) 


850 


49° 35’ 


21° 04’ 


9 


tosie 2 


South (Mountains) 


500 


49° 35’ 


21° 05’ 
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1 


2 


3 


4 


5 


10 


Folusz 


South (Mountains) 


650 


ro 

o 

•''t 


21° 24’ 


11 


Krynica 


South (Mountains) 


700 


49° 21’ 


20° 58’ 


12 


Zakopane 


South (Mountains) 


700 


49° 16’ 


19° 57’ 


13 


Gorce Mount. 


South (Mountains) 


900 


49° 32’ 


20° 05’ 


14 


Swiqtokrzyskie Mount. 


Central (Mountains) 


500 


50° 48’ 


20° 48’ 


15 


Narol 


Central (Uplands) 


345 


50° 20’ 


23° 20’ 


16 


Staszow 


Central (Uplands) 


200 


50° 30’ 


21° 20’ 


17 


Lagow 


Central (Mountains) 


450 


50° 50’ 


21° 10’ 


18 


Brzeziny 


Central (Lowlands) 


200 


51° 50’ 


19° 40’ 


19 


Wielkopolski PN 


Central (Lowlands) 


140 


52° 40’ 


17° 32’ 


20 


Kartuzy 


North (Lowlands) 


200 


54° 15’ 


18° 10’ 


21 


Mlynary 


North (Lowlands) 


120 


54° 00’ 


20° 00’ 


22 


Swierczyna 


North (Lowlands) 


180 


53° 25’ 


16° 15’ 


23 


Gryfino 


North (Lowlands) 


95 


53° 20’ 


14° 45’ 


24 


Swiebodzin 


Central (Lowlands) 


170 


52° 20’ 


15° 20’ 


25 


Sniezka 


South (Mountains) 


755 


50° 50’ 


15° 40’ 


26 


L^dek 


South (Mountains) 


965 


50° 15’ 


16° 50’ 



Tab. 2. Genetic diversity of Polish beech populations 



No. 


Region 


Average 
number of 
individuals 


Average number 
of alleles 
per locus 


Average effective 
number of alleles 
per locus 


Percentage of 
polymorphic 
loci 


Average 

heterozygosity 


Fixation 

index 


Observed 


Expected 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


Cisna 


49.8 


2.1 


1.5 


83.3 


0.250 


0.249 


-0.004 


2 


Moczame 


59.7 


2.5 


1.5 


83.3 


0.241 


0.255 


0.055 


3 


Caryna 


49.7 


2.4 


1.5 


91.7 


0.249 


0.278 


0.104 


4 


Rozsypaniec 


49.4 


2.2 


1.5 


91.7 


0.248 


0.254 


0.024 


5 


Nad Sanem 


57.5 


2.3 


1.5 


83.3 


0.265 


0.269 


0.015 


6 


Hulskie 


58.2 


2.3 


1.5 


83.3 


0.242 


0.250 


0.032 


7 


Wielka Rawka 


50.9 


2.3 


1.4 


83.3 


0.223 


0.244 


0.086 


8 


Losie 1 


49.3 


2.3 


1.5 


83.3 


0.241 


0.262 


0.080 


9 


Losie 2 


45.5 


2.1 


1.4 


91.7 


0.230 


0.229 


-0.004 


10 


Folusz 


45.8 


2.3 


1.5 


91.7 


0.290 


0.308 


0.058 


11 


Krynica 


48.2 


2.2 


1.5 


83.3 


0.266 


0.280 


0.050 


12 


Zakopane 


49.9 


2.4 


1.5 


91.7 


0.268 


0.274 


0.022 


13 


Gorce 


59.7 


2.3 


1.4 


91.7 


0.218 


0.231 


0.056 


14 


Swi^tokrzyskie 


47.4 


2.1 


1.5 


83.3 


0.271 


0.268 


-0.011 


15 


Narol 


48.4 


2.2 


1.5 


83.3 


0.259 


0.272 


0.048 


16 


Staszow 


49.8 


2.0 


1.5 


75.0 


0.242 


0.255 


0.051 


17 


Lagow 


49.7 


2.1 


1.5 


75.0 


0.263 


0.271 


0.030 


18 


Brzeziny 


49.4 


2.1 


1.5 


91.7 


0.233 


0.264 


0.117 
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1 


2 


3 


4 


5 


6 


7 


8 


9 


19 


Wielkopolski PN 


50.5 


2.1 


1.4 


75.0 


0.223 


0.237 


0.059 


20 


Kartuzy 


48.6 


2.0 


1.4 


75.0 


0.224 


0.236 


0.051 


21 


Mlynary 


49.3 


1.9 


1.5 


66.7 


0.242 


0.269 


0.100 


22 


Swierczyna 


49.8 


1.9 


1.4 


75.0 


0.236 


0.249 


0.052 


23 


Gryfino 


46.8 


2.3 


1.5 


83.3 


0.246 


0.268 


0.082 


24 


Swiebodzin 


48.6 


2.0 


1.5 


66.7 


0.222 


0.258 


0.140 


25 


Sniezka 


48.8 


2.3 


1.5 


91.7 


0.254 


0.262 


0.031 


26 


L^dek 


47.8 


2.2 


1.5 


83.3 


0.270 


0.277 


0.025 



Tab. 3. Variation of analysed genetic parameters within different European regions 



No 


Region 


Average number 
of individuals 


Average number 
of alleles per locus 


Percentage of 
polymorphic loci 


Average heterozygosity 


Fixation 

index 


Observed 


Expected 


27 


West Slovakia 


657.3 


2.8 


91.7 


0.269 


0.280 


0.039 


28 


East Slovakia 


424.3 


2.8 


91.7 


0.258 


0.263 


0.019 


29 


Romania 


615.1 


3.3 


100.0 


0.246 


0.263 


0.065 


30 


Ukraine 


123.8 


2.4 


91.7 


0.251 


0.285 


0.119 


31 


Czech + Moravia 


272.1 


2.5 


91.7 


0.254 


0.271 


0.063 


32 


Poland 


1310.3 


3.1 


100.0 


0.246 


0.262 


0.061 



Tab. 4. Characteristics of investigated alleles in Poland and neighboring regions in Europe 



Alleles 


Region 


Observed 


Not observed 


Lap-IA 


Sudety Mts,West Pomerania 


Carpathians, Uplands of Central Poland 


Got-2C 


Carpathians (Bieszczady), Romania, Gryfino 


Sudety, Uplands of Central Poland 


Got-2D 


Carpathians (Bieszczady), West Slovakia 


Other regions 


Idh-C 


Carpathians 


Other regions 


Mdh-IB 


Swi^tokrzyskie Mts., Romania, 


Other regions 


Mdh-ID 


Sudety Mts. 


Other regions 


Mdh-2B 


Other regions 


Sudety 


Mdh-2D 


Romania, Czech + Moravia 


Poland, Slovakia 


Mdh-2E 


Ukraine 


Other regions 


Mdh-3C 


Carpathians, Romania 


Other regions 


Mdh-3D 


Romania 


Other regions 


Pgi-IA 


Carpathians, West Slovakia, Czech, Ukraine 


Sudety, Uplands of Central Poland, Pomerania 


Pgi-ID 


Carpathians (Bieszczady), West Slovakia 


Other regions 


Pgm-A 


All other regions, Carpathians (Bieszczady) 


Other regions of Poland 


Skdh-A 


Carpathians 


Other regions 


Skdh-D 


Slovakia 


Other regions 
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Caryna 


3 


Nad Sanem 


5 


Kartuzy 


20 


Swierczyna 


22 


Cisna 


1 


Moczarne 


2 


Wielkopolski 


19 


Sniezka 


25 


Czech & Moravia 


31 


Rozsypaniec 


4 


Wielka Rawka 


7 


Hulskie 


6 


tosie 1 


8 


Swi^tokrzyskie 


14 


Narol 


15 


Staszow 


16 


Gryfino 


23 


tagow 


17 


Mtynary 


21 


East Slovakia 


27 


West Slovakia 


28 


Gorce 


13 


L^dek 


26 


Krynica 


11 


Zakopane 


12 


Brzeziny 


18 


tosie 2 


9 


Folusz 


10 


Ukraine 


30 


Swiebodzin 


24 


Romania 


29 




0 

Fig. 1 . Dendrograme based on Nei (1972) genetic distances 

cies and observed heterozygosities (Tab. 2) were cal- 
culated based on diploid genotypes. Gene diversities 
(expected heterozygosities) were calculated according 
to Nei (1974) using BIOSYS 1 program (Swofford and 
Selander 1981). Effective number of alleles per locus 
(Crow and Kimura 1970) and fixation index (Wright 
1965) were also calculated to estimate deviation from 
panmictic equilibrium in each population. To assess 
genetic differentiation patterns, dendrogram based on 
Nei (1972) genetic distances was constructed (Fig. 1). 
The results were compared to the genetic diversity 
parameters of beech populations in the neighbouring 
European countries (Tab. 3, Tab. 4) as published by 
Gomory et al. (1995, 1998, 2003). 



0,011 Genetic distance 



Results 

Polish populations were characterized by a slight de- 
crease in the average number of alleles per locus and 
percentage of polymorphic loci towards the north- 
eastern natural range of beech (Tab. 2). The highest 
average numbers of alleles per locus were found in 
the mountain populations from the Sudety and the 
Carpathians Mountains (south Poland): Moczarne 
(2), Caryna (3), and Zakopane (12). Low levels of al- 
lelic richness were characteristic for the lowlands (in 
north, west and central Poland). The highest percent- 
age of polymorphic loci was found in southern Poland 
for populations Caryna (3), Rozsypaniec (4), Losie 2 
(9) in the Carpathians and Sniezka (25) in the Sudety 
Mountains. Unexpectedly high polymorphism was ob- 
served in Brzeziny (18), located on the border of the 
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natural beech limit in Poland, while the lowest value 
was found in Mlynary (21), and in the probably planted 
population Swiebodzin (24). 

The highest values of observed and expected het- 
erozygosity were recorded in the Carpathian popu- 
lations Folusz (10), and Gorce (13). The differences 
between expected and observed heterozygosity were 
the highest in the populations Brzeziny (18), and 
Swiebodzin (24). 

The fixation index was found to be negative in 
beech populations from Cisna (1), Losie 2 (9), and the 
Swietokrzyskie Mountains (14), which can be explained 
by the excess of heterozygotes in the populations and 
some disturbances of typical mating system for bal- 
anced populations in Hardy - Weinberg equilibrium. 
A comparison of genetic parameters from different Eu- 
ropean countries is shown in Tab. 3. The level of genetic 
variation and differentiation estimated in Poland was 
comparable to other regions of Central Europe. Differ- 
ences in the average number of alleles per locus, and av- 
erage heterozygosity were negligible between regions. 
Analysis of the frequency distribution of different al- 
leles allowed for the division of regions characterized 
by the presence or absence of certain alleles (Tab. 4). 
This distribution of alleles might reflect the migration 
paths during the Holocene. 

Discussion 

The genetic differentiation of Polish beech populations 
is comparable to those in other neighbouring Europe- 
an countries including Slovakia, the Czech Republic, 
Ukraine, and even Romania, which was reported previ- 
ously (Paule et al. 1995). 

Generally, alleles can be divided due to the range 
of alleles of large-scale example Mdh-1 as well as genes 
encoding peroxidases: Px-1 and Px-2. Some alleles were 
locally or regionally specific, namely Lap-IA (Sudety 
Mountains, Pomerania), Got-2C (Bieszczady Moun- 
tains, Gryfino, Romania), Mdh-IB (Swiqtokrzyskie 
Mountains, Romania), Mdh-3D (Romania), Mdh-3D 
(Ukraine), Pgi-ID (Bieszczady Mountains, West Slova- 
kia), Skdh-D (Slovakia). 

The dendrogramme (Fig. 1) constructed based on 
the Nei’s (1972) genetic distances does not allow dis- 
crimination of regional beech populations. The dendro- 



gramme indicates small genetic differentiation between 
populations geographically located far from each other 
such as Caryna (3) and Nad Sanem (5), compared with 
Kartuzy (20) and Mlynary (21). 

The highest genetic variation between populations 
was found in the East Carpathians. There is a slight de- 
crease in genetic polymorphism to the north of Poland, 
which can be explained by postglacial migration paths 
from south to north (Szafer 1935; Raiska -Jasiewiczowa 
1983; Hazier et al. 1997 as well as Gomory et al. 1999, 
and Magri et al. 2006). 

It turns out that it is difficult to find specific mark- 
ers for beech provenances that allow for discrimina- 
tion between geographical regions, due to the observed 
small genetic differentiation. This may be explained by 
the restricted occurrence of this species in small refu- 
gia during the glacial period with subsequent expansion 
into the currently occupied range during a short time 
(about 5.000 years). High genetic variability and the 
presence of alleles Lap-1 A and MdhlD not occurring in 
the Carpathians, but only in Sudety populations might 
confirm this opinion, which may support the observed 
expansion of beech to Polish territories from one refu- 
gium (Magri et al. 2006). 

Beech from the north could slightly increase its 
range and connect to that migrating from the south. 
Raiska- Jasiewiczowa (1983) indicates that about 1.000 
BP it was noticed the significantly large participation of 
beech in the Swiqtokrzyskie Mountains in comparison 
to surrounding regions. The combination of beech refu- 
gium from the south and migrating from north along 
the lake may have occurred just in the region, which 
may explain the rare occurrence of allele Mdh-IB in the 
region. 

Conclusions 

- Polish beech populations are characterized by high 
genetic diversity, similar to other European beech 
populations. 

- Towards the north of the natural range limit, the 
beech populations are poorer in allelic diversity and 
are genetically less well differentiated than those in 
southern regions. 

- The highest genetic differentiation was found in 
the East Carpathians (Bieszczady Mountains). The 
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analysis of paths of beech migration after the glacial 
period from one refugium for Central Europe could 
be an explanation of this phenomenon. 
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